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We study the photon generation in a transmission line oscillator coupled to a driven qubit in the 
presence of a dissipative electromagnetic environment. It has been demonstrated previously that a 
population inversion in the qubit may lead to a lasing state of the oscillator. Here we show that the 
circuit can also exhibit the effect of "lasing without inversion". This is possible since the coupling 
to the dissipative environment enhances photon emission as compared to absorption, similar to the 
recoil effect which was predicted for atomic systems. While the recoil effect is very weak, and so 
far elusive, the effect described here should be observable with present circuits. We analyze the 
requirements for the system parameters and environment. 
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Introduction. The basic "circuit quantum electro- 
dynamics" (cQED) system consists of a superconduct- 
ing qubit coupled to a transmission line oscillator. The 
former replaces the atom, the latter the radiation field 
of the traditional quantum electrodynamics setupi. In 
recent experiments^ a cQED version of the single-atom 
maser was realized. For that purpose a superconducting 
single-electron transistor is coupled to a transmission- 
line resonator with frequencies in the microwave regime. 
In this particular case a population inversion of the two 
charge states, which are in resonance with the oscillator, 
is achieved by pumping via a third stated"—. The same 
setup can also be used to create non-classical photon 
states in the resonator— or, when coupled to a mechanical 
oscillator, it is one of the prime candidates to create and 
observe non-classical states in macroscopic objects^ii^. 

The exchange of energy quanta between the atom and 
the cavity, or in the present case qubit and resonator, in- 
cluding the escape of photons from the cavity, is governed 
by the balance 
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Here P„ is the probability to find n photons in the cavity, 
while P-|- /I are the occupation probabilities of the atomic 
levels. The rates of stimulated photon emission and ab- 
sorption are given by T^j^ and respectively, while 
K is the decay rate of photons in the oscillator with n 
thermal photons. In thermal equilibrium the occupation 
probability decreases as a function of the photon num- 
ber, Pn+i/Pn = n/{n + 1). In contrast, the lasing state 
is characterized by a sharp peak of P„ at nonzero values 
of n, and hence P„+i > Pn below the peak. In the opti- 
cal domain the absorption and emission rates are roughly 
equal. Hence the condition that P„ should grow with n 
can only be achieved if a population inversion, P-|- > Pj^, 
is created. This is done in a pump process which typically 
involves a third level. Examples are the early microwave 
lasers based on ammonia moleculesii as well as the recent 



cQED lasers with charge qubits, where the pumping is 
achieved by quasiparticle tunneling processes via a third 
state with an extra quasiparticle^. 

As lasing technology progressed it became clear that 
population inversion is not really needed. One of the 
earliest schemes of lasing without inversion (LWI) is the 
pumping of Rabi sideband a^^d^ ; its cQED analog was 
demonstrated and discussed in Refs. [ij and [H. How- 
ever, this process is not strictly LWI, since closer inspec- 
tion reveals that it is based on a population inversion in 
the dressed states basis^. LWI without hidden inversion 
can only be achieved by breaking the symmetry of pho- 
ton emission and absorption^^. The only experimental 
demonstrations of this symmetry breaking were accom- 
plished by using several external fields, which change the 
absorption and emission profile of the bare ato m^^d^ . 

In this paper we will employ a fundamentally differ- 
ent method that is more similar to the idea of using 
the shift of the emission peak caused by the recoil ef- 
fect on the atom^°. In this case the atom always ab- 
sorbs energy for both emission and absorption. This ef- 
fect has never been verified experimentally in conven- 
tional lasing setups, since it requires large frequencies 
and/or low temperatures such that the recoil energy ex- 
ceeds the temperature.— With cQED systems the situ- 
ation appears more promising due to two differences to 
quantum optics. One reason is that there are many in- 
trinsic sources of non-classical noise suppressing energy 
emission as compared to absorption. The other is the 
strong coupling between the qubit and the resonator, 
which makes it possible to create a lasing situation even 
in the presence of strong noise. 

In the following we will introduce the system and the 
methods needed to describe the shift of the emission 
peak. Our methods reproduce the standard lasing 
results^^ but, in addition, describe the effect of strong 
noise as long as the photon number is not to large. 
We will then demonstrate that it is possible to use 
a noisy qubit to create a non-thermal photon state 
in the oscillator and that the oscillator field exhibits 
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the linewidth narrowing characteristic for the lasing 
state. We describe the requirements for the system 
and environment needed to observe inversionless lasing. 
Finahy, we wiU describe a specific experimental setup 
based on a superconducting circuit. 

The System: The basic mechanism which leads to 
lasing without inversion is the shift of the photon emis- 
sion line due to the coupling to a dissipative environ- 
ment. This effect may play a role in every setup that can 
be described by the extension of the Jaynes— Cummings 
Hamiltonian {h = 1), 

H = wofl' a + g{a)(T- + aa+) + ]^/\Ea^ + ^X{t)az ■ (2) 

Here a is the annihilation operator for the oscillator mode 
with frequency wq, the Pauli matrices ai act on the states 
I t) and I I) , the qubit energy splitting and the strength of 
the coupling between oscillator and qubit are denoted by 
AE and g, respectively. We also account for fluctuations 
X{t) of the energy splitting of the qubit. They may origi- 
nate from the electrodynamic environment and, together 
with strong coupling, make inversionless lasing possible. 
The noise X{t) is characterized by its correlator, 
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where T is the base temperature. To be specific we as- 
sume that the fluctuations result from an Ohmic envi- 
ronment; taking into account a capacitive coupling to 
the system it leads to a spectral function of the form 



J(aj) = ec w X 
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parametrized by a coupling strength ec and cut-off fre- 
quency UJR. 

Here we consider situations where the noise is strong. 
Therefore we proceed using the polaron transformation 

U = exp[— zcr^ y_^X{t')dt'], which yields 

H = Ho + Hf, 
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In the following we will assume that the Hamiltonian 
part Hg leads to incoherent transitions, which we ana- 
lyze in the spirit of the so-called P(£')— TheoryS-^ by an 
appropriate Liouville term in the Master equation. This 
somewhat unusual treatment of the coupling between 
oscillator and qubit differs from the approaches to lasing 
followed, e.g., by Refs. S [ll [H, and [13 where the 
coupling is treated as part of the coherent time evolu- 
tion. However, as we will demostrate below our master 
equation approach is valid in a broad parameter regime. 



and it reproduces results for the average photon num- 
ber and phase coherence time which are typical for lasing. 

Master Equation. Our system is described by the 
Master equation 



P = i[Ho, p] + [/^diss + -Cg + £pump] P ■ 
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It contains three Liouville operators, corresponding to an 
incoherent pump process, -CpumpP, photon emission and 
absorption, £gP, and dissipation in the oscillator, CdissP- 
Incoherent pumping is described by 



-CpumpP^ rup(2o-+p(7_ - [(T_0-+,p] + ) 

+rdown (2(7- pcr+ - [o-+(T-,p] + ) 



(7) 



where [, ]+ denotes an anti-commutator. Without cou- 
pling to the oscillator, g = 0, the occupation probabil- 
ities of the atom, P-j- = {f \p\ and similar Pj,, would 
be given by Pf/Pj, = Tup/rdown- That means, as long 
as Fdown ^ Tup the pumping does not produce a popula- 
tion inversion. For later use we introduce the parameter 
Tpump — Tup + Tdown and define the population inversion 

COefficent Dq = (Fup - Fdown )/rpuinp- 

The crucial effect to be considered here is the emission 
and absorption of photons, i.e., the transitions between 
the states | t)!*^) and | i)|n -I- 1) induced by the Hamil- 
tonian Hg. To calculate the rates we expand the time 
evolution of the density matrix up to second order in 
and consider the level broadening caused by the pump. 
In this case the Liouvillian is given by 

F+ 

£gP ~ — ^ (2a^a^ paa^ — [aa+a^a^ , p]+) (8) 
F~ 

_l — |h (2aa+pa-a^ - [a-a^aa+, p]+) , 

with photon emission and absorption rates Fpj^ = 
<?^'S'ph(±(5w). They depend on the spectral function 



^ph(c^)- / Cph(t)e-r— l*l/2e-*df, (9) 
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at the frequency given by the detuning Suj = AE — uq. 
We assume that a fairly strong pump is applied to the 
qubit, as it was the case in the experiments of Ref. j^. In 
this case the lowest order expansion in Hg converges for 
all combinations of ec and ujn., as long as g\/n <^ Fpump. 
Including the level broadening caused by the pump is cru- 
cial to reproduce the standard lasing results in the limit of 
weak noise ec — >■ 0, ujr — >■ oo. One should note that even 
for g^/n > Fpu^np our master equation is still valid for 
£(7 = in the stationary limit and for g^n <C \/kBTec^R 
at all times. 
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Finally, dissipation in the oscillator is described by the 
standard Lindblad operator^"', 

AissP = -^(^ + 1) {2apa^ ~ [a'^a,p]+) (10) 

+ |n(2aVa- [aa^p]+) , 

with dissipation strength determined by the decay rate 
K and the temperature entering the thermal number of 
photons n — (exp(a;o/fc_BT) — 1)^^. 

The master equation ^ can be solved by numerically 
diagonalizing it in the eigenbasis of Hq. On the other 
hand, there are two standard methods to solve the 
Master equation analytically^^. One method is to 
consider the equation of motion of operator averages and 
the second method is to adiabatically eliminate qubit 
states and derive an effective equation of motion for the 
oscillator photon number states, pnn' = {n\T r quhit p\n' ) , 
where the trace is taken over the qubit eigenstates | f) 
and I 4,). With the combination of these two methods 
we can derive analytical results for all quantities we are 
interested in. To find a closed set of equations for the 
operator averages we have to make the semi-classical 
approximation (a^aaz) ~ {a^a){(Jz)- This is a good 
approximation in the stationary limit and allows us to 
derive an analytical expression for the average photon 
number in the system (n) = (at a). To describe the time 
evolution of the system we use the adiabatic elimination 
of the qubit degrees of freedom. This method accurately 
describes the time-evolution of the off-diagonal matrix 
elements even for ^ T^^. 

Lasing Without Inversion. To demonstrate that 
inversionless lasing is possible in cQED we show the av- 
erage photon number in the oscillator for three illustra- 
tive cases (see fig. [T]). We always choose the detuning 
between oscillator and qubit such that the photon emis- 
sion rate Fpj^ is maximized. The cases that are shown 
correspond to the standard lasing situation with a noise- 
less qubit, 6(7 —> 0, LJR — ?> cxD, and to situations with 
strong noise. Without noise the emission and absorption 
rates are equal and we can only achieve heating of the 
oscillator. For strong noise we tune the system such that 
there is additional energy available if a photon is created, 
5u) = AE — w > 0. The additional energy is given to the 
environment. This means that for the qubit to absorb 
a photon the environment has to deliver the energy Suj. 
Since we consider an environment at low temperatures, 
we get strong suppression of photon absorption and hence 
there is a strong imbalance between the emission and the 
absorption rate, Fpj^ F",^. So despite the fact that 
no population inversion is created in the qubit a large 
number of photons is excited in the oscillator. For the 
noiseless case the number of photons only starts to rise 
strongly as we achive inversion, Dq > 0. 

The optimal condition for inversionless lasing is 
reached if the coupling to noise is stronger than the level 
broadening caused by the pump ec ^ Fpump and we 
maximize the asymmetry due to small temperatures by 



0.02 



0.01 



pump 




-1 1 

Su/ec 




FIG. 1. a) The photon emission and absorption rate, Fp^^ 
and Fpjj, as a function of the frequency detuning Slj. With- 
out noise, lur — >■ oo, the two rates are equal (black line). 
Dashed, red and blue lines show the emission and absorption 
rate, respectively, for Fpump/ec = 0.0325 and ujR/ec = 0.5. 
Full, red and blue lines show the rates for Fpump/ec = 0.0325 
and ujn/ec = 0.05. b) The average photon number (n) as 
function of the inversion coefficient Dq. The black line is 
the result for a qubit without noise, ua — >■ oo, at Su = 0. 
The red lines are the results under the influence of noise, for 
u)R = 0.5 at 5uj/ec = 0.05 (dashed line) and for ujr = 0.05 
at Siu/ec = 1 (full line). The frequency detuning has been 
chosen to maximize the photon number. We used the pa- 
rameters Fpump/ec = 0.0325, ksT/ec = 0.05, g/ec = 0.01, 
n/ec = 8.125 X 10"^ 



choosing ec ^ ksT > lor. For each photon that is 
emitted or absorbed the energy ecr^ is given to the en- 
vironment. Therefore the maximum for photon emission 
is at Sbj = ec where photon absorption is strongly sup- 
presed. The rates Q obey detailed balance only for the 
system in equilibrium, Fpu^p — 0- However, if a pump 
is applied the extra broadening has to be taken into ac- 
count. At the maximum of the photon emission rate we 
can approximate the rates by 



ph 



(11) 



From these rates we can formulate the threshold condi- 
tion defined in eq. ([1]) for n ^ 1 as 

(12) 
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Using the explicit form of the photon emission and ab- 
sorption rate (|lip we will now summarize the two cru- 
cial prerequisites for inversionless lasing: strong noise at 
small temperatures and strong coupling. The first con- 
dition is needed to create an imbalance between photon 
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FIG. 2. a) The phase correlation decay rate func- 
tion of the frequency detuning Suu. The values of the cut- 
off frequency are ivn/ec ~ 0.05 (red), iVR/ec = 0.5 (green), 
^r/^c = 2 (blue), b) The average photon number (n) as a 
function of the frequency detuning 5a; with color coding as 
in a). Full lines display the results of a numerical solution of 
the master equation, dashed lines are the approximate ana- 
lytical solution p5p . The system parameters are Fpump/ec ~ 
0.0325, ksT/ec = 0.05, g/ec = 0.01, w/ec = 8.125 x 10-^ 
Do = -0.25 . 



emission and absorption. Comparing the absorption and 
emission rate given by eq. (jlip we find 
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The second condition, strong coupling, can be formulated 
by comparing photon emission with the oscillator decay 
rate. This yields 
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One should note that the condition of strong coupling 
requires only a comparison of coupling to noise strength. 
We still want the system to be within the limit of the 
rotating wave approximation, <C w, in contrast to the 
so called ultra-strong coupling regime recently observed 
in Ref. HI. 

For the conditions of strong noise (|13|) and strong 



coupling (jl4p we additionally assumed y/ksTec 3> ujjj. 
However this is not strictly necessary as is shown in fig. [2] 
Here we show the average photon number obtained from 
a numerical calculation of the photon emission and ab- 
sorption rates. Even for larger cut-off frequencies we can 
create a lasing state in the oscillator. However, the pho- 
ton number gets maximal for a small cut-off frequency. 

For unequal photon creation and annihilation rates and 
assuming that all conditions formulated in eqns. (1121) . 
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to be 



14]) apply we find the average photon number 
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This results confirms our previous discussion. To get a 
high photon number we need a low decay rate k in the 
oscillator and FupFpj^ — rdownr~jj > 0. We also see the 
standard lasing condition emerge. For equal photon cre- 
ation and annihilation rates, Tpj^ = we would need 
Tup > Tdown to get a large photon number. If all orders 
of K are taken into account our result for the average 
photon number exactly reproduces the standard lasing 
results in the limit ec — >■ 0,Wi{ — ^ oo. 

One of the characteristic properties of a laser is the 
phase coherence, which is expressed in the small decay 
rate of the phase correlator C(t) = (a(i)a(O)). In the 
Markovian limit the correlator decays exponentially with 
the rate kl, C{t) cx e~'^'^*. To find the correlator we use 
the regression theorem 



(a^(t)a(O)) =Tr[aVaW] 



(16) 



where Pait) is a density matrix which follows the time 
evolution given by the master equation (|6]) with the initial 
condition /9o(0) = apst- After adiabatically removing the 
time-evolution of the qubit, the problem is reduced to the 
decay of the matrix elements pnn+i- For large photon 
number the decay rate is given by 
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Even for Fpj^ 7^ Fpj^ the decay rate is inversely propor- 
tional to the photon number as it is the standard case for 
a laser. In the limit of weak noise, ec ~^ 0; ~^ 00, this 
result reproduces the standard lasing results, as long as 
we stay within the general limit of validity of our master 
equation. For a discussion of the phase correlation decay 
rate in the limit of strong coupling and weak noise, see 
e.g., Ref. 0. 

In fig. [5] we show the phase correlation decay rate kl 
and the corresponding average photon number (n) as a 
function of the frequency detuning Suj. The numerical 
and analytical results for the average photon number are 
in good agreement; for the decay rate the comparison 
is not shown but there is also qualitative agreement. 
Results are shown for different values of the cut-off 
frequency w^j, and one can see that the enhancement of 
photon number is the stronger, the smaller is w^j. On 
the other hand, even for larger values of ujr lasing is 
possible with a strongly reduced phase correlation decay 
rate. 

Outlook. A specific system where the conditions for 
lasing without inversion can be realized is the single 
artificial- atom laser investigated by Astafiev et al.^. This 
system consists of a charge qubit coupled to an oscilla- 
tor. An applied transport voltage induces quasiparticle 
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tunneling via a third state which can be used to create 
a population inversion in the qubit, which has led to the 
observed lasing. With the same setup it would also be 
possible to only create an enhanced population of the ex- 
cited qubit level without inversion. The noise spectrum 
needed for inversionless lasing can be created by coupling 
the charge qubit to an external resistor (see e.g. Ref. [26h . 
To reach a low cut-off frequency it is necessary to use a 
resistor with a large resitance, as it has been demon- 



strated, e.g., in Ref. H^. The strong coupling needed to 
satisfy the requirement is standard in most cQED 
experimentsii^i^. 
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